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A key actor in the conventional theory of superconductivity is the induced interaction between 
electrons mediated by the exchange of virtual collective fluctuations, originally phonons. Other 
collective modes that can play the same role, especially spin-fluctuations, have been widely dis¬ 
cussed in the context of high-temperature and heavy Fermion superconductors. The strength of 
such collective fluctuations is measured by the associated susceptibility. Here we use differential 
elastoresistance measurements on five optimally doped Fe-based superconductors to reveal that a 
diverging nematic susceptibility appears to be a generic feature in the optimal doping regime of 
these materials. The observation motivates consideration of the effects of nematic fluctuations on 
the superconducting pairing interaction in this family of compounds, and possibly beyond. 

PACS numbers: 74.70.Xa, 72.15.-v 


A growing body of evidence suggests the possibil¬ 
ity of an intimate connection between electronic ne¬ 
matic phases [T] and high-temperature superconductiv¬ 
ity. However, it is currently unclear to what extent 
there is any causal relationship between nematic fluc¬ 
tuations and superconductivity. Strongly anisotropic 
electronic phases have been found in the underdoped 
regime of both cuprate EM j and Fe-based EHE5 high- 
temperature superconductors. For underdoped Fe-based 
systems, recent measurements of the elastoresisistance 
[IMS, Raman spectroscopy p^j - fTH] . and elastic mod¬ 
uli [T9j [20! for the representative electron-doped system 
Ba(Fei_ x Co x ) 2 As 2 reveal a divergence of the electronic 
nematic susceptibility upon approach to the tetragonal- 
to-orthorhombic structural phase transition, definitively 
establishing that the phase transition is driven by elec¬ 
tronic correlation. For the cuprates, recent x-ray diffrac¬ 
tion j2H - [24j and NMR [25 measurements have revealed 
evidence for short-range charge density wave order in 
“underdoped” crystals. Although details of the charge 
ordered state(s) are still being established, these initial 
observations have at least motivated discussion of a pos¬ 
sible “vestigial” nematic order [26]. Perhaps significantly, 
in the phase diagrams of both families of compounds, op¬ 
timal doping is located close to putative quantum critical 
points P71 - E5] which potentially have a nematic character 

EH HUGS]. 

From a theoretical perspective, recent treatments in¬ 
dicate that nematic quantum criticality, (i.e. quantum 
critical fluctuations caused by proximity to a nematic 
quantum critical point) can provide an enhancement of 
the existing pairing interaction. In particular, a pure ne¬ 
matic phase does not break the translational symmetry of 
the original crystal lattice; consequently the q =0 nematic 
fluctuations enhance T c in all symmetry channels l-'l 1 II.‘{.'S| . 


It is therefore of considerable interest to empirically es¬ 
tablish whether nematic fluctuations are a characteristic 
feature of optimally doped high temperature supercon¬ 
ductors, as well as to probe the extent to which these 
nematic fluctuations show intrinsically quantum behav¬ 
ior. In the present paper, we show that this is the case 
for Fe-pnictide and chalcogenide superconductors by con¬ 
sidering the representative materials Ba(Fei_ a; Co a; ) 2 As 2 
and Ba(Fei_ a; Ni a; ) 2 As 2 (i.e. electron-doped “122”), 
Bai_ a; K a ,Fe 2 As 2 (hole-doped), BaFe 2 (Asi_ a ,P a; )2 (isova- 
lent substitution) and FeTei_a;Se x (“11”). Furthermore, 
we find that the nematic susceptibility obeys a simple 
Curie-Weiss power law for all five optimally doped Fe- 
based superconductors over a wide temperature range. 
For the electron and hole-doped 122 pnictides a sub- 
Curie-Weiss deviation was observed at low temperatures, 
which we tentatively attribute to an enhanced sensitivity 
to disorder in a quantum critical regime. 

Nematic order couples linearly to anisotropic strain 
of the same symmetry. Consequently, the nematic sus¬ 
ceptibility of a material can be measured by considering 
the electronic anisotropy that is induced by anisotropic 
in-plane strain. In the regime of infinitesimal strains, 
all forms of electronic anisotropy are linearly propor¬ 
tional. Hence, the rate of change of resistivity anisotropy 
with respect to anisotropic strain, defined in the limit of 
vanishing strain, is linearly proportional to the nematic 
susceptibility[M]. The proportionality constant depends 
on microscopic physics, but away from any quantum crit¬ 
ical point does not contain any singular behavior. Con¬ 
sequently, the induced resistivity anisotropy reveals the 
essential divergence of the nematic susceptibility upon 
approach to a thermally-driven nematic phase transition 
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For a tetragonal material, nematic order has either d xy 
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FIG. 1: Temperature dependence of the f? 2 9 elastore- 
sistance of BaFe2As2. The data follow a Curie-Weiss be¬ 
havior, which is the anticipated mean-field temperature de¬ 
pendence of the nematic susceptibility of a material approach¬ 
ing a thermally driven nematic phase transition I13H15I . Up¬ 
per panel shows —2m^&, proportional to the nematic suscep¬ 
tibility Xn(b 2 )• Black line shows the Curie-Weiss fit. The 
quality of fit can be better appreciated by considering the 
inverse susceptibility, — (2mee — 2mg 6 )^ 1 which is perfectly 
linear (left axis of lower panel; fit shown by red line), and the 
Curie constant — (2m^ — 2mg 6 ) _1 (T — T*) (right axis of lower 
panel), which is independent of temperature. The Weiss tem¬ 
perature obtained from the Curie-Weiss fit, which gives the 
bare mean field nematic critical temperature, yields a value T* 
= 109 ±1.4 K. Coupling to the lattice renormalizes the critical 
temperature, leading to a nematic-structural phase transition 
at T s = 134 K (vertical gray line). Inset shows a photograph 
of a square crystal glued on a PZT piezoelectric stack for dif¬ 
ferential elastoresistance measurements using Montgomery’s 
geometry. Four electrical contacts were made at the corners 
and a strain gauge was glued on the top surface [36]. 


symmetry (i.e. the B 2g irreducible representation of the 
D^h point group, corresponding to nematic order ori¬ 
ented along a nearest-neighbor Fe-Fe bond - the [110] or 
[110] crystal axes), or d x 2_ y 2 symmetry (Bi gi correspond¬ 
ing to the [100] or [010] crystal axes). Anisotropic strains 
with appropriate symmetry are then £6 = (e^y + £yx)/2 
and ei — £2 = £ ra — Cyy respectively. The strain-induced 
changes in resistivity can be described using the dimen¬ 
sionless elastoresistivity tensor, m.y m- 
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(A p/p) i = ^m i je j (1) 

j =i 

where 1 = xx, 2 = yy etc. Consequently, the correspond¬ 
ing components of the nematic susceptibility tensor are 
given by: 


Xn(b 2 „) = c x 2ro 66 (2) 

XN(B lg ) =c' x (mn - m 12 ) (3) 

where c and c' are proportionality constants, which de¬ 
pend on microscopic physics [51], 

We measure elastoresistivity by applying an in-situ 
tunable anisotropic strain using a piezoelectric PZT 
stack. A square plate sample (typical dimension 
750x750x20 pm) is glued on the side wall of the PZT 
stack using a commercial two part epoxy. The PZT stack 
deforms when a voltage is applied and hence strains the 
sample glued on top of it [35]. The amount of strain can 
be measured by a strain gauge glued either on the back¬ 
side of the PZT stack or on the top surface of a larger 
sample(the latter case enabling a full determination of 
the strain transmission[36l). The in-plane resistivity ten¬ 
sor of the sample is measured via the Montgomery tech¬ 
nique, with electrical contacts made at the four corners of 
the square sample [36] . Representative data taken using 
this new technique are shown in Figure 1 for the specific 
case of BaFe 2 As 2 - As has been previously demonstrated 
mm, the data can be fit very well by a Curie-Weiss 
temperature dependence: 


2m 66 = 2mg 6 + A/[a(T-r*)]. (4) 

In Figure 2 we show B 2g elastoresistance data 
for a range of optimally doped materials, in¬ 
cluding BaFe 2 (Asi_ x P a; )2 (isovalently substituted), 
Ba(Fei_ x Ni a; ) 2 As 2 (electron-doped), Bai_ a; K 2 ,Fe 2 As 2 
(hole-doped) and FeTei-^Se^,. In all cases, 2?Ti66 rises 
strongly with decreasing temperature, with comparably 
large values for each compound. The observation of such 
an effect for this wide variety of ways of doping, including 
the case of the iron chalcogenide FeTeo. 6 Seo. 4 , is highly 
suggestive that a divergence of the nematic susceptibil¬ 
ity in the B 2g channel is a generic feature of optimally 
doped Fe-based superconductors. This is our main re¬ 
sult. Regardless of microscopic models, from a purely 
empirical perspective it is apparent that the optimally- 
doped superconductor is born out of an electronic state 
that is characterized by strongly fluctuating orientational 
(nematic) order in this specific symmetry channel. This 
result is especially notable for FeTe 0 . 6 Se 0 .4 given that the 
orientation of both the magnetic ordering wave vector 
and the in-plane component of the structural distortion 
of undoped FeTe is 45 degrees away from the orientation 
of those in the iron arsenide materials. 

The 77166 coefficient of optimally doped 
BaFe 2 (Aso. 68 Po. 32)2 (Fig. 2A) follows a perfect 

Curie-Weiss temperature dependence from T = 250AT 
(the highest temperature where strain can be effectively 
transmitted by the epoxy) down to T c (below which 
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FIG. 2: Divergence of the f? 2 9 elastoresistance 2meo of several different families of optimally doped Fe pnic- 
tide and chalcogenide superconductors. (A) optimally doped BaFe 2 (Aso.rPo. 3)2 (isovalent substitution), (B) optimally 
doped Ba(Feo.955Nio.o45)2As2 (electron doped), (C) optimally doped Bao.5sKo.42Fe2As2 (hole doped), and (D) optimally doped 
FeTeo. 58 Seo. 42 . Insets indicate the dopant site (red) in the respective unit cells of each material. Comparable magnitude induced 
resistivity anisotropies are observed for each case j36l . Upper panels show |2m66|, whereas lower panels show |(2m,66 — 2mg 6 )| -1 
(left axes of lower panels, blue symbols) and |(2m-66 — 2m6 6 )|(T — T*) (right axes of lower panels, black curves). Black(upper 
panels) and red(low panels) lines shows fits to Curie-Weiss behavior of m 66 and |(2m66 — 2 to 6 6 )|~ 1 respectively. Grey horizontal 
lines (low panels) shows the average values of |(2m66 — 2m66)|(T — T*) in the fitting temperature range. Regions of deviation 
from Curie-Weiss behavior in (B) and (C) are indicated by gray shaded regions. For (A) and (B), 2m66 is negative. For (C) 
and (D), 2mee is positive. Fit parameters are listed in I36ll44l . 
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FIG. 3: Variation of the f? 2 9 elastoresistance of Ba(Fei_ a ,Co a: ) 2 As 2 for four representative compositions. (A,B) 
underdoped compositions Ba(Feo.97sCoo.o25)2As2 and Ba(Fe 0 .953Coo.o47)2As2; (C) optimal-doped Ba(Fe 0 .g3Co 0 .07)2As2; and 
(D) overdoped Ba(Feo. 92 Coo.os) 2 As 2 . The heavily underdoped composition is very well described by a Curie-Weiss temperature 
dependence over the entire temperature range(black lines in upper panels). For the compositions near optimally-doping, 2mse 
can be well fit by a Curie-Weiss T-dependence at high temperatures. Below a characteristic temperature scale(different for 
each compositions, indicated by shade gray region), a strong downward deviation from Curie Weiss behavior is observed, also 
seen in the inverse susceptibility (upward curvature) and in |(2m66 — 2m6 6 )|(T — T*) oc Xn(T — T*) (strong downturn), which 
are shown in the lower panels. The deviation from Curie-Weiss behavior is the strongest at the optimal doping, and diminishes 
on either side of the phase diagram. Curie-Weiss fit parameters for each composition are listed in (361 . 
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temperature the resistance is zero), with a Weiss tem¬ 
perature (T*) close to zero. For FeTeo.eSeo. 4 , the W66 
coefficient can be perfectly fitted with Cure-Weiss down 
to T c but shows a significant deviation for temperatures 
greater than 100K, possibly related to the loss of quasi 
particle coherence due to its extremely small Fermi 
energy as observed in photoemission spectroscopy and 
transport]!?, [38]. Intriguingly, for the electron and 
hole doped 122 pnictides, a downward deviation from 
Curie-Weiss behavior is observed below a characteristic 
temperature that is different for each of the materials. 
The quality of fit to the Curie-Weiss functional form for 
BaFe 2 (Aso. 68 Po. 32)2 in comparison to all of the other 
optimally doped compositions studied can be readily ap¬ 
preciated if the data are plotted on log axes (Fig. Sll). 
Only the data for the P-substituted system can be fit by 
a single power law [|(2 to 66 — 2mg 6 )| ~ (T — T*) -7 ] over 
the entire temperature range, yielding 7 = 0.985 ± 0.005 
and T* = 11.7 ±3.1. 

A slightly overdoped BaFe2(Aso.64Po.36)2 sample has 
also been measured; the elastoresistivity coefficient m^Q 
of this sample can also be well fitted by Curie-Weiss tem¬ 
perature dependence over the entire temperature range, 
with a negative Weiss temperature T* = — 11.5±2.3(Fig. 
S14). The small value of the Weiss temperature observed 
in BaFe 2 (Asi_ a ,P a; )2 near optimally doping, and the fact 
that it crosses zero as doping increases, motivates con¬ 
sideration of the nematic susceptibility in the context of 
a quantum critical point. An Ising nematic phase transi¬ 
tion in an insulator would generally be expected to have 
dynamical exponent z = 1 , so the effective dimensional¬ 
ity of the system would be d + z = 3+1. As the materials 
in question are metallic, the situation is more compli¬ 
cated [39]; the Hertz-Millis paradigm is expected to break 
down in the case of electronic nematic order in metallic 
systems, for which Landau damping associated with the 
gapless Fermions leads to a larger value of the dynamical 
exponent, 2 = 3. Although the analysis of this problem 
is not straightforward, both analytic [40, and numeri¬ 
cal studies m indicate that the nematic susceptibility 
at criticality should diverge in proportion to 1/T (up to 
possible logarithmic corrections), consistent with the ob¬ 
served behavior of the B 2g elastoresistance of optimally 
doped BaFe 2 (Asi_ rc P a ;) 2 . 

To gain insight into the physical origin of the low 
temperature downward deviation from Curie-Weiss be¬ 
havior in the electron and hole doped pnictides, we 
consider the evolution of the B 2g elastoresistance as 
a function of composition for the specific case of Co¬ 
doped BaFe 2 As 2 - Figure 3 A-D show the progression 
of 2m66 as a function of Co doping for several repre¬ 
sentative compositions (more were measured, see [36]), 
from underdoped (A,B), through optimal doping (C), 
to the overdoped regime (D). For the heavily under¬ 
doped compositions(x = 0.025), 2?7i66 can be well de¬ 
scribed by mean field Curie-Weiss T-dependence down to 


T s . It has been previously established that this behav¬ 
ior is essentially independent of disorder, at least com¬ 
paring undoped BaFe 2 As 2 with different Residual Resis¬ 
tance Ratios (RRR), and Co and Ni substituted crys¬ 
tals with the same TV m- The Weiss temperature T* 
extracted from high-temperature fits to Curie-Weiss be¬ 
havior crosses zero close to optimal doping (Fig. 4). 
However, for lightly underdoped Ba(Feo. 553 Coo.o 47 ) 2 As 2 , 
a downward deviation from mean-field behavior at low 
temperatures begins to be noticeable, and becomes more 
pronounced for optimally doped Ba(Feo. 93 Coo.o 7 ) 2 As 2 ; 
similar deviations are observed in optimally Ni and K 
doped BaFe 2 As 2 - A similar effect has been observed 
in recent measurements of the sheer modulus via three- 
point bending experiments [90] , The deviation from 
Curie-Weiss behavior diminishes again as the doping is 
further increased: for the overdoped composition (Fig. 
3D) the data can be fit to Curie-Weiss to a lower tem¬ 
perature than for optimal doping, and the magnitude 
of the deviation below ~ 45 K is smaller than for opti¬ 
mal doping. 277766 °f four additional Co dopings has also 
been measured, showing a similar non-monotonic doping 
dependence of the deviation from Curie-Weiss behavior. 
The fact that the effect is maximal near optimal doping 
where T* ~ 0, suggests that it is associated with prox¬ 
imity to the putative nematic quantum phase transition. 

P-substituted BaFe 2 As 2 is the “cleanest” of all of the 
known 122 Fe-pnictide families, evidenced by the fact 
that quantum oscillations can be observed across the 
phase diagram [42] [43]. The deviation from Curie-Weiss 
behavior for optimally doped compositions of all other 
dopants in BaFe 2 As 2 suggests that disorder plays an im¬ 
portant role in the quantum critical regime. All forms 
of quenched disorder produce locally anisotropic effec¬ 
tive strains, which thus couple to the orientation of the 
nematic order; this is “random-field” disorder [26]. Anal¬ 
ysis of the random-field Ising model yields several gener- 
ically expected effects of random field disorder, includ¬ 
ing suppression of the nematic susceptibility below mean 
field expectations for a clean system, and, for the case of 
a quantum phase transition, the enhanced sensitivity of 
quantum critical phenomena to disorder. pISj. 

The temperature dependence of the nematic suscep¬ 
tibility can also be extracted from measurements of the 
elastic moduli :|5D|. The two measurements (elastoresis¬ 
tance and elastic moduli) are in broad agreement, for 
example in terms of the Curie-Weiss T-dependence of 
Xn for underdoped compositions of Ba(Fei_ a; Co a ;) 2 As 2 , 
and also in terms of the deviation from Curie-Weiss be¬ 
havior near optimal doping. However, there is an im¬ 
portant distinction in terms of the relative magnitude of 
the measured quantities as a function of doping. In par¬ 
ticular, the normalized lattice softening [cg 6 — cqq]/c%q 
extracted in m monotonically decreases in magnitude 
and extent in temperature as a function of x. In con¬ 
trast, the quantity | (277X66 — 2m.g 6 ) I = c Xn initially in- 
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FIG. 4: Phase diagram of Ba(Fei_ I! Co a ;)2 Ast, showing 
the variation of 2niQg in the x — T plane (color scale). 

Squares, circles and triangles indicate T s , Tjv, and T c re¬ 
spectively. Stars indicate the bare mean-field nematic critical 
temperatures, T* , extracted from the Curie-Weiss fits of 2mse 
above the temperature at which disorder effects suppress the 
nematic susceptibility (see text; light grey symbols are used 
for the cases for which deviations from Curie Weiss behavior 
are observed at low temperatures). As has been previously de¬ 
termined via longitudinal elastoresistance measurements [13]. 
but established here by the full i? 2 S differential elastoresis¬ 
tance, the Weiss temperature, T* goes through zero as a func¬ 
tion of x close to optimal doping. Color scale shows the mag¬ 
nitude of —2rri66, which peaks between lightly underdoped to 
optimal doped compositions. 

creases with x, peaking for lightly underdoped composi¬ 
tions x — 0.05. The apparent enhancement of the elas¬ 
toresistance 277166 over the softening of the elastic mod¬ 
ulus for compositions near optimal doping is potentially 
related to renormalization of the quasiparticle effective 
mass in the quantum critical regime, as has been seen 
in P-substituted BaFe 2 As 2 [29] ■ Such a mass renormal¬ 
ization is an expected consequence of nematic quantum 
critical fluctuations [39]. Significantly, it is precisely these 
low-energy quasiparticles that are also involved in the 
eventual superconductivity. 

The temperature dependence of the nematic suscep¬ 
tibility can also be extracted from measurements of the 
elastic moduli m- The two measurements (elastoresis¬ 


tance and elastic moduli) reveal broad agreement, for ex¬ 
ample in terms of the Curie-Weiss T-dependence of \n 
for underdoped compositions of Ba(Fei_ 2; Co a: ) 2 As 2 , and 
also in terms of the deviation from Curie-Weiss behav¬ 
ior near optimal doping. However, there is an impor¬ 
tant distinction in terms of the relative magnitude of the 
measured quantities as a function of doping. In partic¬ 
ular, the normalized lattice softening [cg 6 — c 66]/c66 ex- 
tracted by Boehmer et al. [2D] monotonically decreases 
in magnitude and extent in temperature as a function 
of x. However, as can be seen from inspection of Fig¬ 
ure 3, the quantity |(2?ti,66 ~ 2rrie 6 ) I = c Xn initially in¬ 
creases with x, peaking for optimally doped composi¬ 
tions. The apparent enhancement of the elastoresistance 
2m66 over the softening of the elastic modulus for com¬ 
positions near optimal doping is potentially related to 
renormalization of the quasiparticle effective mass in the 
quantum critical regime, as has been seen for instance 
in P-substituted BaFe 2 As 2 [29]. Such a mass renormal¬ 
ization is an expected consequence of nematic quantum 
critical fluctuations [59]. Significantly, it is precisely the 
low-energy quasiparticles that determine the elastoresis¬ 
tance that are also involved in the eventual superconduc¬ 
tivity. 
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Materials and Methods 

Elastoresistivity coefficients 

Elastoresistivity relates changes of resistivity, p , of a 
solid conductor to strains that it experiences. We define 
the change in resistivity with respect to a strain e as 

mm- 


Ap a p = p a p(e) - Pa/3(e = 0) (5) 


And the normalized resistivity change as mm- 


(Ap/p) a p 



( 6 ) 


Since resistivity and strain are both 2nd rank tensors, 
the elastoresistivity, M, is a 4th rank tensor defined by 
the following expression [45] : 


(A-Pjp) otp ^ ^ 

7 S 

a, 13,7,6 = x, y, z 


(7) 

( 8 ) 


In the absence of a magnetic field, the Onsager reciprocal 
relation holds for the electrical conductivity, and the as¬ 
sociated symmetry properties of the resistivity and strain 
tensors allow us to more conveniently express them as six 
component vectors, such that 

6 

(A p/p)k = m ki e i (9) 

/=! 


where 1 = xx, 2 = yy,3 = zz, 4 = yz,5 = zx, 6 = xy. 
Using this Voigt notation, the elastoresistivity is repre¬ 
sented by a pseudo-second rank tensor. For a tetragonal 
crystal, the point group symmetry allows us to further 
simplify the elasotresistivity tensor: 


/ mu mi 2 mi 3 0 0 0 \ 

mi2 mu nii3 0 0 0 

m 3 1 m 3 i m 33 0 0 0 

0 0 0 m.44 0 0 

0 0 0 0 m 44 0 

y 0 0 0 0 0 m-66 / 


( 10 ) 


Here the six compenent vectors describing the resistivity 
change and the strain are defined with respect to the a, 
6, and c crystallographic axes which are aligned with the 
x , y 1 and z Cartesian axes. The elastoresistivity coeffi¬ 
cients will be written with respect to the crystallographic 
axes throughout this text. In the following sections, we 
describe how the elastoresistivity coefficients (mn — 77112 ) 
and 7B-66 are related to the nematic susceptibility for the 
Big and -B 2g symmetry channels, and how these coeffi¬ 
cients can be measured. These ideas were first described 
in ref [2] ■ 





Elastoresistivity coefficient and nematic 
susceptibility 


For a C 4 to C 2 phase transition the nematic order pa¬ 
rameter is a Ising variable, and is proportional to the 
anisotropy of any equilibrium electronic properties U|. 
The dimensionless resistivity anisotropy 


N = Pxx Pvv (11) 

\Pxx + Pyy) 

does not purely measure the equilibrium electronic prop¬ 
erties, but is determined by the anisotropy in the elec¬ 
tronic structure and the anisotropy in the scattering rate; 
the latter being a dynamical property of the system. Nev¬ 
ertheless, one can still express the nematic order param¬ 
eter $ as a function of the resistivity anisotropy N and 
take a Taylor expansion for small value of N, for which 
the leading term is linear. Therefore, for small value of 
N, which applies to our measurements (which are all in 
the linear regime), it can be related to the nematic or¬ 
der parameter via a proportionality constant N = c^>. 
Strain with the same symmetry (i.e. e xx — e yy ) couples 
linearly to dq and hence can be considered as a symme¬ 
try breaking field conjugate to the electron nematic order 
parameter. Therefore we can define the nematic suscep¬ 
tibility xn as the quantity that relates the induced order 
parameter and symmetry breaking strain in the region 
of linear response. For a tetragonal crystal there are two 
possible symmetry channels in which nematic order can 
develop via a continuous phase transition. The two sym¬ 
metry channels correspond to Bi 9 and B 2g irreducible 
representation of a tetragonal point group, which couples 
linearly to orthorhombic distortions oriented along [100]t 
and [llOj-r direction in the tetragonal lattice. Hence for a 
crystal oriented with [100]t along the x direction (define 
by the PZT stack orientation): 

XN(Blg) = d ^Blg = c ( mn _ m 12 ) (12) 

d{e xx -e yy ) 

Similarly, for a crystal oriented with [llOj-r along the x 
direction: 


d^B2g , , 

XN(B2g) — -Ie—0 — C TO 6 6 


(13) 


Therefore the elastoresistivity coefficient is related to 
the bare nematic susceptibility via the same proportion¬ 
ality constant c. The proportionality constant c depends 
on the microscopic details of the electronic structure, the 
precise definition given to the nematic order parameter, 
and, in some cases, the nature of the disorder scatter¬ 
ing, and in general could be temperature and materials 
dependent. However, two experimental observations of 
iron based superconductors suggest that c is a constant 
within a wide range of temperatures and independent 
to a certain degree of disorder[15]. First, the temper¬ 
ature dependence of the elastoresistivity coefficient mee 


of the parent compounds BaFe 2 As 2 can be perfectly fit¬ 
ted by a Curie-Weiss function (Fig. 1 of main paper), 
which is exactly the same functional form of \n derived 
by a Ginzburg-Landau theory. Second, the elastoresis¬ 
tivity coefficient mee of BaFe 2 As 2 has been measured for 
crystals with different residual resistance ratios (RRR), 
which is a measure of crystalline disorder [15] ■ It was 
found that mee is independent of the RRR of the sample. 
A similar result is found for underdoped compositions of 
Co- and Ni-doped BaFe 2 As 2 which have the same Neel 
temperature [To] . The perfect fitting to Curie-Weiss in 
these conditions strongly suggest that c remains a con¬ 
stant over a wide range of temperatures and different 
degrees of disorder. 


Elastoresistance Measurements using a modified 
Montgomery method 

The elastoresistivity coefficients (mu — 777 - 12 ) and mee 
can be determined from measurements of p xx and p yy as 
a function of anisotropic strain e xx — e yy [TU [TS] . Previ¬ 
ously, we have demonstrated that anisotropic strain could 
be achieved by gluing thin crystals to the side surface of a 
PZT piezoelectric stack [46] . The resistivity p xx and p yy 
were measured separately by gluing two long bar samples 
along the poling and the transverse directions of PZT 
stack. In this paper we introduce an improved technique 
by measuring p xx and p yy on a single square sample via a 
modified Montgomery method. The advantage of using a 
single square sample is twofold: by using a single sample 
we remove the uncertainties due to unequal strain expe¬ 
rienced by two crystals; the square shape also eliminates 
possible errors due to geometric effects, ensuring that the 
strain transmission is equal along the x and y directions. 
In the remainder of this section, we describe the details 
of this new technique. 

Single crystals were grown via standard self flux tech¬ 
niques, as described elsewhere m ■ Square samples were 
cut from larger crystals such that the tetragonal c-axis 
was perpendicular to the plane of the thin squares, and 
the a-axis aligned either along the sides of the square 
(suitable for probing the B lff channel) or at 45 degrees 
(for B 2 c,). Anisotropic strain was achieved by gluing these 
thin crystals to the side surface of a PZT piezoelectric 
stack]46j. with the sides of the squares aligned along the 
principle x and y axes of the stack. By applying a posi¬ 
tive voltage bias through the voltage leads, the PZT stack 
causes the sample to experience an anisotropic strain 
by expanding along the poling direction and contracting 
along the transverse direction. We define the poling di¬ 
rection of the PZT stack as the y-direction and the trans¬ 
verse direction as ^-direction. The resistances along the x 
and y-directions of a single sample were measured simul¬ 
taneously while the voltage applied to the PZT stack was 
varied (typically for three full hysteresis loops between - 
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50 V and +150 V) (Fig. SI). A photograph and illustra¬ 
tion of the configuration of a typical experiment probing 
the elastoresistance of such a sample is shown in Figures 

51 and S2. Current is sourced at a given frequency across 
the top two contacts A and B of the sample, while the 
resistance along the ^-direction (R x ) is measured by the 
voltage drop across the bottom two contacts C and D. 
Simultaneously a second current at a separate, distinct 
frequency is sourced through the leftmost two contacts, 
C and D, while the resistance along the y-direction {R y ) 
is measured across the rightmost contacts, D and B. 

The resistivity along the x and y-directions p X xiPyy 
can be calculated from the measured resistances R x , R y 
using a modified Montgomery method if the sample 
side lengths (L x and L y ) and thickness (L z ) are known 
(Fig. S 3 ) HU 05 ]. The square anisotropic crystal can 
be mapped to an isotropic equivalent sample with side 
lengths L' x and L' y using the expression: 


Bl 

L' 


; ta £; + l/ ( i ta £> ,+4 


(14) 


where R x and R y are the measured resistances along the 
x and y- axes respectively [181 05 ] . In the thin sample 
limit (———r < 0.5) this ratio can then be used to cal- 

(L*L y )2 

culate the resistivities along the PZT axes according to 
the relations 0H10S3 : 





L 7 L t . . . L 


7 tL' 


L n 




L 'v 


(15) 

(16) 


Neglecting geometric effects (ie small changes in L x 
and L y of the strained sample, which are discussed in 
the next section), p xx and p yy can be obtained for each 
applied strain at each temperature. Now we can express 
the experimentally measured resistivity change for a thin 
square sample experiencing B\ g strain (e.g. it is glued so 
the square edges and the tetragonal a-axis of the crystal 
are aligned to the axis of the PZT stack) in terms of the 
elastoresistivity tensor by using the following expression: 


(A p/p) xx = e xx mu + e yy m 12 + e zz m 13 (17) 

= e xx (rnu - vpmu ~ vsmi 3 ) (18) 

(A p/p) vv = e xx (m 12 - z'pmii - ^mi 3 ) (19) 


to obtain: 

(A p xx - A p yy )/p = (e xx - e yy )(mn - mi 2 ) (20) 
The normalization factor, p = ^/p xx p yy , is calculated 
for each sweep using p xx and p yy values taken at zero 
anisotropic strain. If the dynamic range of the PZT is not 
large enough to tune through the zero anisotropic strain 
point the zero strain values are determined by linearly 
extrapolating p xx and p yy to the value at which they are 
equal. Since zero anisotropic strain does not necessarily 
coincide with zero isotropic strain, the B\ g response and 
the normalization factor are determined at a point with 
some residual isotropic A\ g strain. 

Expression S 16 shows that mu — m \ 2 is the appropri¬ 
ate coefficient for the B ly nematic susceptibility. Sim¬ 
ilarly, by cutting the crystal in a square configuration 
with the edges oriented at 45 degrees with respect to the 
tetragonal a-axis (so now it is experiencing B 2g strain 
with respect to the crystal axes), it can be easily shown 
via transformation of the elastoresistivity tensor that the 
difference in the change of resistivities can be expressed 
by the following: 

(A p xx ^Pyy) / P — (Cxx C yy )2?7l66 (21) 

Here x and y still refer to the axes along which normal 
strain is applied and m66 is defined with respect to the 
crystal axes. Since e xx can be measured directly by the 
strain gauge glued on the back of the PZT stacks[50| and 
the Poisson ratio is known ]14j . we can directly extract 
the elastoresistivity coefficients by fitting the linear slope 
of (Ap xx — Ap yy ) / p as a function of strain. An example of 
the linear fitting can be found in Fig. S4. Further details 
and the associated description of the elastoresistance can 
be found in refs. mm 


Geometric factor 

During an elastoresistivity measurement the sample 
undergoes a small physical deformation, therefore in ad¬ 
dition to the change of resistivity that is proportional to 
the elastoresistivity coefficients, there is also a change of 
the sample dimensions (ie a geometric effect that also 
affects the measured resistance). In a standard four ter¬ 
minal resistivity measurement this can be expressed as: 

Ap/p = AA/A — AL/L + AR/R (22) 


where up = —e vv /e xx is the Poisson’s ratio of the PZT 
stack and us = —tzz/^xx is the effective Poisson’s ratio 
of the sample glued to the stack. We take the difference 
of the change of the resistivity of the two measurements 


Where A is the cross sectional area of the current path 
and L is the separation between contacts. For the Mont¬ 
gomery method the equivalent expression is more com¬ 
plicated and is given by: 









10 


Ap/p = A A/A — AL||/L|| 


foil A^-L - R ± AR\\)Coth(^ ^Z 2+ln{ *" ] ) (M%) + y/** 2 + M|fo 2 ) 

27ri?|| R ± ^JAtt 2 + ln(^) 2 
2tt( - fl||Ai*L + i?j.Ai?||(1 + y / 4 

27ri?|| i?j_ y/47r 2 + in(^) 2 


(23) 


where Ly is the length between the measurement leads, 
i?y is the resistance of the sample along the current path, 
and R± is the resistance of the sample perpindicular to 
the current flow. For an isotropic square sample (i?y = 
Rj_ = R and Ai?y = ARj_ = A R) this simplifies to the 
standard four terminal expression. 

In both cases, the geometric terms A L/L and A A/A 
account for the change in the resistance of the crystal due 
solely to changes in physical dimensions and have no tem¬ 
perature dependence. In a simple metal, such as copper, 
the elastoresistivity coefficients are small, and the change 
of resistance is dominated by the geometric factor. How¬ 
ever, in semiconductors or strongly correlated materials 
the elastoresistivity coefficients can be much larger than 
the geometric effect I5U- In the analysis described in the 
main text, the geometric factor is partially responsible 
for the small, temperature independent term . 


Practical considerations: Strain Transmission 

Strain is transmitted to the sample through a stan¬ 
dard two part epoxy[52] used to glue the samples on to 
the sidewall of the PZT stacks. Strain can be dissipated 
both in the epoxy layer and in the sample itself. The 
ability of the stack to transmit the strain to the sample 
therefore depends on the choice of glue and the thickness 
and elastic properties of the sample [53]. Previously it was 
found that the strain is essentially fully transmitted for 
large (i.e. in plane dimensions >2000^m) samples with 
thicknesses below approximately 100/rm [I3j ; measure¬ 
ments presented in this paper were made for samples with 
thicknesses <50/rm. To confirm that the strain was trans¬ 
mitted in our samples (~ 750/im square), we prepared a 
series of test samples: three square BaFe 2 As 2 samples of 
similar thickness with side lengths ~ 300/irn, ~ 750/xm, 
and ~ 3000/mi. Strain gauges [501 were glued both on the 
back side of the stack and on top of the largest test crystal 
(due to the large size of the strain gauge we are not able 
to glue the gauge on top of the smaller crystals). In Fig. 
S5 we show the strain transmission in the largest crystal. 
For temperatures below 255K, the strain transmission is 


> 80% with only a weak temperature dependence. We 
compare the m§Q elastoresistivity coefficient calculated 
from the strain measured from the two (sample and PZT 
mounted) strain gauges in Fig. S6. The two calcula¬ 
tions are in good agreement with each other and have 
Curie-Weiss critical temperatures (T*) values of 109. OK 
± 0.7 K and 107. 7K ± 0.5 K respectively, where uncer¬ 
tainties are standard errors obtained from each individ¬ 
ual Curie-Weiss fit. The weak temperature dependence 
of the strain transmission seen in Fig. S5 is dominated by 
the much stronger Curie-Weiss temperature dependence 
of the elastoresistance and therefore has negligible effect 
on the determination of the elastoresistivity coefficient 
m.66 and the associated fit parameters. Details about the 
Curie-Weiss fit and fit parameters can be found in the 
main text and subsequent sections. 

Though we cannot directly measure the strain in the 
smaller crystals we can compare the 777-66 strain response 
(Fig. S7). We see good agreement in the high tempera¬ 
ture m66 response between the ~ 3000/rm and ~ 750/iin 
samples, indicating a similar strain transmission in the 
smaller ~ 750/rm sample and larger ~ 3000/iin sam¬ 
ple. However, the smallest sample, while looking qual¬ 
itatively similar with a divergent nematic susceptibility, 
has a smaller response. This indicates that the strain is 
not fully transmitted, but we can still compare the tem¬ 
perature dependence of the response. By normalizing 
the 77i66 by its value at 135A' ( mee(l35K)), we see that 
the normalized responses of the two samples are identical 
within uncertainties (Fig. S8). The fitted T* values are 
107.3/t ± 0.9 K and 109.0 K ± 0.7 K for the small and 
large samples respectively. The similarity in the tem¬ 
perature dependence indicates that the small sample be¬ 
haves as if it was simply experiencing a smaller strain, 
which does not effect the extracted nematic critical tem¬ 
perature. While the majority of samples used in this 
paper are in the size range of full strain transmission, 
it is important to note that the phosphorus substituted 
iron arsenide samples are ~ 300/im square and are in 
the regime of smaller strain transmission. While this af¬ 
fects the quantitative estimate of ni66 for these samples, 
it does not affect the temperature dependence. See Table 
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SI for a list of all sample dimensions. Comparison of the 
Weiss temperature extracted from these three measure¬ 
ments of BaFe 2 As 2 samples (109 .OK ± 0.7 1\ for 3000/im 
sample, 116 K ± 2.2 K for 750/im sample, 107.3A' ± 0.9AT 
for 300/im sample) providing an estimate of the sys¬ 
tematic uncertainties that presumably comes from vari¬ 
ation in sample adhesion, which for example could affect 
the amount of residual Ai 9 strain that could potentially 
change the system’s bare nematic critical temperature. 

Comparison to differential longitudinal 
elastoresistance using two crystals 

Additional data were collected using the differential 
longitudinal elastoresistance technique, in which two bar 
like samples were used to obtain p xx and p yy separately, 
following the method we originally described in ref Ca¬ 
in comparison to the newly developed modified Mont¬ 
gomery’s method that only requires single square shape 
sample, the bar shape sample suffers from strain relax¬ 
ation along the direction perpendicular to the current 
if the width of the sample is comparable to its thick¬ 
ness. Nevertheless for materials with a elastoresistance 
dominated by the response in nematic channel, these two 
techniques yield qualitatively similar results. 

To demonstrate this point, data collected by differen¬ 
tial longitudinal technique are shown in Fig. S9 (four dif¬ 
ferent optimally doped iron-based superconductors), and 
Fig. S10 (four dopings of Ba(Fei„ a ,Co a; ) 2 As 2 ). All the 
features presented in Fig. 2 and 3 in the main text includ¬ 
ing diverging nematic susceptibility for optimally doped 
sample, perfect Curie-Weiss temperature dependence for 
P doped BaFe 2 As 2 , and deviation from Curie-Weiss be¬ 
haviors near optimally doping for Ba(Fei_ x Co x ) 2 As 2 are 
reproducible. These agreements further corroborate the 
robustness of our findings. 

Sign and magnitude of m 66 of Bao. 6 Ko. 4 Fe 2 AS 2 

The main subject of our paper is the observation of 
a divergence of the differential elastoresistance that we 
have found for many different optimally doped composi¬ 
tions of Fe-based superconductors. However, beyond this 
main result several other features of the data are of inter¬ 
est in their own right and deserve additional comment. 
In this section we particularly emphasize two aspects of 
the elastoresistance of optimally doped Bao. 6 Ko. 4 Fe 2 As 2 . 
First; whereas our measurements, described in the main 
text, reveal that the sign of the W66 elastoresistivity coef¬ 
ficient is negative for Co, Ni and P substituted BaFe 2 As 2 , 
we find that it is positive for K-doped BaFe 2 As 2 - Second, 
the absolute magnitude of the coefficient is comparable 
for all four materials. Both of these observations warrant 
some additional commentary. 


For Co, Ni and P substituted BaFe 2 As 2 , the negative 
sign of m66 in the tetragonal phase (for temperatures 
above the structural transition temperature, T$) is con¬ 
sistent with the resistivity anisotropy that is observed for 
detwinned crystals in the orthorhombic state for temper¬ 
atures below Ts That is, the resistivity in 

the tetragonal phase decreases as the sample is stretched 
along the current direction. 

For Bai_ x Kj;Fe 2 As 2 , the in-plane resistivity 
anisotropy has also been measured up to 1 = 0.26 
L5S3- It was observed that for x < 0.235, the sign of 
resistivity anisotropy is the same as that of the parent 
compound and Co-doped BaFe 2 As 2 - But for x > 0.235, 
the sign of resistivity anisotropy was found to reverse; 
i.e. the resistivity was found to be higher for the longer 
lattice constant direction. Our observation of a positive 
value for TO66 for optimal-doped Bai_ x K x Fe 2 As 2 which 
has a potassium concentration x = 0.4, confirms this 
dramatic change in sign of the resistive anisotropy for 
sufficient hole doping. Significantly, our measurement 
is made in the tetragonal state, and therefore is not 
complicated by issues associated with the reconstruction 
of the Fermi surface in the antiferromagnetic state, or 
by concerns associated with detwinning crystals in the 
orthorhombic state. 

Despite the agreement in terms of the sign of the re¬ 
sistive anisotropy for sufficiently doped Bai_ x K x Fe 2 As 2 , 
the relative magnitude of the anisotropy of detwinned 
underdoped compositions is quite different to the strain- 
induced elastoresistance near optimal doping. The mag¬ 
nitude of the resistivity anisotropy of detwinned crystals 
of underdoped K-substituted BaFe 2 As 2 was found to be 
rather small [56j . This observation has been interpretted 
by some authors as evidence supporting the notion that 
the resistivity anisotropy derives from anisotropic scat¬ 
tering associated with lattice defects, which acquire an 
anisotropy due to the nematic fluctuations [571. Within 
such a scenario, defects in the FeAs planes are believed 
to have a stronger effect, and hence would naturally yield 
a larger resistivity ansiotropy, than defects that lie away 
from the FeAs plane (as would be the case for K substi¬ 
tution), potentially accounting for the significantly lower 
resistivity anisotropy observed for detwinned underdoped 
Bai_ x K x Fe 2 As 2 relative to Ba(Fei_ x Co x ) 2 As 2 . Our 
measurements, however, reveal a comparable magnitude 
(within a factor of two) for mge of optimally doped Co¬ 
substituted and K-substituted BaFe 2 As 2 , in contrast to 
the much larger difference in the resistivity anisotropy 
in the orthorhombic state for underdoped compositions. 
Clearly K substitution per se does not necessarily imply 
a small resistive anisotropy. 

There are three possible explanations for the small 
anisotropy that was previously measured in the or¬ 
thorhombic state for K-substituted BaFe 2 As 2 - Firstly, 
it is possible that the resistivity anisotropy has a small 
magnitude simply because it changes sign as a function of 
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composition, and therefore must go through zero for some 
composition. Secondly, it is possible that the FS recon¬ 
struction in the antiferromagnetic state affects the resis¬ 
tivity anisotropy differently for electron and hole doped 
cases. Finally, it is also possible that the samples used for 
the study of underdoped compositions were not fully de- 
twinned, although this seems less likely given the careful 
characterization of the twin populations by the authors 
[55] , At present, it is not clear which of these possibili¬ 
ties is responsible for the small resistivity anisotropy of 
underdoped samples of Bai_ a ,K a; Fe 2 As 2 in the antiferro¬ 
magnetic orthorhombic state. 

Bearing this in mind, P-substitution gives another use¬ 
ful insight to the origin of the resistive anisotropy. Be¬ 
yond the residual resistivity, another good measure of the 
disorder affecting transport in the Fe plane is the pres¬ 
ence or absence of quantum oscillations. P-substituted 
samples exhibit quantum oscillations across the phase di¬ 
agram (at least on the overdoped side[42l [43],) indicat¬ 
ing a much lower elastic scattering rate than the other 
compounds we have discussed so far, for which quan¬ 
tum oscillations are not observed. Significantly, niee of 
the P-substituted samples is comparably large to that of 
Co, Ni and K doping, despite the lower elastic scattering 
rate, implying that the physical origin of the dominant 
contribution to the resistive anisotropy is unlikely to be 
associated with scattering from anisotropic defects, con¬ 
trary to suggestions advanced in ref. 153, even though 
such effects are necessarily present. 

The observations described above serve to emphasize 
an important point; that elastoresistance measurements 
for T > T s are the most relevant data to compare with 
theories of the resistivity anisotropy in the nematic phase, 
rather than comparison with transport anisotropies in 
the low T Neel and/or orthorhombic state measured in 
samples detwinned by uni-axial stress. There are several 
reasons, but primarily because (a) there are no concerns 
about degree of detwinning (in the tetragonal phase the 
material is homogeneous, and homogeneously strained); 
(b) there are no issues with back-folding of bands and FS 
reconstruction; and (c) we have precise knowledge of the 
strain experienced by the material (which is not known 
for detwinned samples, for which the amount of strain de¬ 
pends on the temperature dependent elastic modulus). In 
other words, in selecting transport data to compare with 
theoretical predictions for resistive anisotropy, elastore¬ 
sistance measurements provide the more quantitatively 
controlled experimental perspective. 

Zero strain resistivity versus temperature of the 
optimal doped Fe-based superconductors 

For reference, we show in Figure Sll, the temperature- 
dependence of the resistivity for the same materials for 
which elastoresistance data are shown in the main text. 


The data were taken for free-standing (unstrained, and 
unattached) crystals. Data are plotted as R/R(300K) to 
eliminate uncertainty in geometric factors. 

Elastoresistivity coefficients of Fei+aTei _ :r Se„, 

The ”11” family of iron chacogenides (Fei+aTei-^Sea,) 
provide an interesting example of the richness of the 
phase diagram of iron based superconductors. At the 
tellurium end of the phase diagram, Fei+^Te possesses 
a different antiferromagnetic ground state; the antifer¬ 
romagnetic order is a bicollinear order with an order¬ 
ing wave-vector that is 45° to that in iron pnictides)58|. 
Consequently the spontaneous strain developed in the 
symmetry breaking phase is e xx — e yy (Bi g ) rather than 
e xy (B 2g ). On the other hand, at the selenium end of 
the phase diagram, FeSe only has a structural transi¬ 
tion that develops a spontaneous B 2 g strain e xy without 
any long range magnetic order [F5]. The optimal doped 
FeTeo. 6 Seo .4 we measured in this work has a T c = 13A' 
and there is no sign of phase transitions other than the 
superconducting transition. 

We have measured mu — m 12 (B lg ) elastoresis¬ 
tivity coefficient of FeTeo.6Seo.4, Bao.6Ko.4Fe2As2 and 
Ba(Feo. 93 Coo.o 7 ) 2 As 2 - These measurements were taken 
using a differential elastoresistance technique which in¬ 
volves simultaneously measuring the resistivity of two or¬ 
thogonal bar-like samples. Details of this technique can 
be found in ref. m- The results are plotted in Fig. S12 
together with the mee (B 2 g ) coefficient presented in Fig. 
S9. It can be clearly seen that the nematic suscptibil- 
ity diverges only in the B 2 g channel rather than the Bi g 
channel. 


Further details about the Curie-Weiss fit parameters 

As described in the main text, mee data were fit to a 
Curie-Weiss temperature dependence: 

2m 66 = 2m° 6 + A/[o(T - T*)] (24) 

where T* gives the bare mean-field nematic critical 
temperature (the temperature at which nematic order 
would occur if there were no coupling to the crystal lat¬ 
tice Hi TO), and A/a is the Curie constant. 

The fitting is performed by varying T* and 2mg 6 un¬ 
til the standard deviation of (2mee — 2mg 6 )(T — T*) is 
minimized, in which 2m66 and T are experimentally mea¬ 
sured data. This fitting method is justified by the fact 
that (2m66 — 2m,g 6 )(T — T*) = A/a is expected to be a 
constant as a function of temperature. 

The determination of fitting temperature range is not 
a trivial task, there are systematic deviations from Curie- 
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Weiss behavior in both low temperatures (due to the ef¬ 
fect of disorder near quantum critical point) and high 
temperatures (due to the vanishing nematic fluctua¬ 
tions). In order to unbiasedly determined the high and 
low temperature cutoff for the fitting, we have developed 
a standardized procedure: A fitting is first performed 
over the widest possible temperature range in which the 
strain is still effectively transmitted, i.e. from T$ (for 
underdoped samples) or Tq (for optimal and overdoped 
samples) to 250A'. The T* and 2rng 6 determined by 
this fitting is then used to calculate the deviation of 
(2?7i66 — 2mg 6 )(T — T*) from its mean value as a func¬ 
tion of temperature. An example for optimally doped 
Ba(Fei_ a; Co a ;) 2 As 2 sample is shown in the top panel 
of Fig. S13. There are systematic deviations in both 
low and high temperature ranges, which is typically ob¬ 
served in all the data sets. The next step is to lower 
the higher temperature cutoff to the deflection point in¬ 
dicated by the arrow in Fig. S13, then perform the fit¬ 
ting in the new temperature range. The deviation of 
(2?n 66 — 2?7ig 6 )(T — T*) is calculated again, as shown in 
the second panel of Fig. S12. In the case of optimally 
doped Ba(Fei_ a; Co; C ) 2 As 2 , the systematic deviations has 
been reduced yet they are still observable. Therefore the 
next step is to increase the lower temperature cutoff un¬ 
til the systematic deviations is smaller than the intrinsic 
noise level, as shown in the third panel of Fig. S13. Fi¬ 
nally, we raise the upper temperature cutoff again until 
the systematic deviation appears to be greater than the 
intrinsic noise level (4th panel in Fig. S13). 

For heavily underdoped Ba(Fei_ ;c Co a ;) 2 As 2 , and for 
optimally doped BaFe 2 (Aso. 7 Po. 3)2 and Fei+^Tei^Sea,, 
all the systematic deviations can be eliminated at the 
second stage of the standardized procedure, i.e. there 
is no need to raise the lower temperature cutoff, and 
?7Zgg follows a perfect Curie-Weiss temperature depen¬ 
dence all the way down to T s or T c . For optimally doped 
Co, Ni and K-substituted BaFe 2 As 2 , the data can be 
described by Curie-Weiss behavior only above a temper¬ 
ature that is higher than T c . The range of temperatures 
over which this results in a reasonable description of the 
temperature-dependence of in gg is shown by the linear fit 
(red line) to the inverse susceptibility |2 ?ti 66 — 2wg 6 |~ 1 
in the lower panels of Figures 2 and 3 in the main text. 
Four additional dopings of Ba(Fei_ a; Co :c ) 2 As 2 were mea¬ 
sured and the data were used to construct the color map 
phase diagram in Fig. 4. Fig. S14 shows the raw data of 
all eight dopings of Ba(Fei_ 2 ,Co ;c ) 2 As 2 and their CW fit¬ 
tings. Associated fit parameters of all the data presented 
in the main text and supplementary materials are listed 
in Table SI. 

The quality of fit to Curie-Weiss behavior for optimally 
substituted BaFe 2 (Aso. 7 Po. 3)2 can also be appreciated 
when the data are plotted on a log-log scale (Fig. S15). 
As can be seen, only the optimally P-substituted data 
(red) can be fitted linearly over the whole temperature 


range, with a slope 7 = 0.985 ± 0.005 (solid black line). 
For all other optimally-doped samples, the data exhibit 
a deviation from Curie-Weiss behavior. 

In addition, we show below the elastoresistance coef¬ 
ficient 2mg6 of a P-substituted sample which is slightly 
overdoped with respect to the sample shown in the main 
paper. As shown in Fig. S16, for this slightly over¬ 
doped composition the elastoresistance still follows a per¬ 
fect Curie-Weiss temperature dependence, but now the 
Weiss temperature is -11.5 K. For the approximately op¬ 
timally doped composition shown in the main text, the 
Weiss temperature is + 11.7K. Clearly the Weiss tem¬ 
perature crosses zero, indicative of a nematic quantum 
critical point. 

Theory of the Effect of Weak Disorder on the 
Nematic Susceptibility 

Here we analyze the properties of an electron nematic 
that derive primarily from general symmetry consider¬ 
ations. In particular, we discuss how the presence of 
quenched disorder affects the behavior of the nematic 
susceptibility upon approach to a nematic thermal or 
quantum critical point. Admittedly, many aspects of the 
measurements we have reported cannot be understood 
in this way as they depend on material-specific details 
and require an explicit understanding of the processes 
that contribute to the transport anisotropy, specifics of 
the particular type of quenched disorder involved, and 
of course an understanding of the underlying thermody¬ 
namic mechanism of nematicity. Nevertheless, this sim¬ 
ple analysis serves to offer intuition for the effects of dis¬ 
order in nematic systems. It also reveals effects that have 
a remarkable similarity to the observed deviations from 
mean field Curie-Weiss behaviors of the elastoresistance 
for disordered optimally doped compounds presented in 
the main text. 

We start by defining the local nematic order cf>(r ), 
which is even under inversion but odd under rotation by 
7t/2 or reflection through the appropriate mirror plane. 
The uniform thermal average of (f) is the nematic order 
parameter, N = fi” 1 ]TA(0(r)), where f2 = JW is pro¬ 
portional to the volume. In computing the susceptibility, 
we will consider the effect of a symmetry breaking field, 
h , which is conjugate to N , and which thus represents 
the appropriate element of the strain tensor - we define 
the coupling so that positive h favors positive N. Fi¬ 
nally, the quenched disorder is represented by a random 
field b(f) which couples in the same way as h to the ne¬ 
matic order parameter, but is assumed to have average 
strength zero and to be short-range correlated, b(r) = 0 
and b(f)b(f r ) = a 2 5r, r', so a is the relevant measure of 
the strength of the disorder. 

Weak disorder: If the disorder is weak and the sys¬ 
tem is not too close to criticality, the effect of disorder 
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on the nematic susceptibility can be computed perturba- 
tively in the disorder strength, which yields a result of 
the form 

N = X (T,a)h-Xs(T,cr)h 3 (25) 


X" 1 (T, a) = x” 1 (T) + a 2 T(T) + ..., (26) 

where x(T,0) = x(T), Xs(T,0) = X 3 (T), and T{T) are 
all related to equilibrium properties of the disorder free 
system (i.e. a = 0, for which formally exact expressions 
can be readily obtained in terms of various correlation 
functions of (f >), and ... refers to higher order terms in 
powers of h and <r. 

We have defined signs in these equations such that un¬ 
der typical circumstances Xi X 3 j and r are all positive 
quantities. To obtain explicit expressions, we can con¬ 
sider the problem in the context of a Landau-Ginzberg 
free energy of the form 

f [ 4 >] =\ 0(r) G_1 ( f ' - 

r> " (27) 

r r 

which we treat as a mean-field functional which we solve 
by minimizing F with respect to <f>. In the disordered 
phase (where N —> 0 as h —► 0), this gives rise to the 
expressions 

X(T) = gi, 

r(T) = 3ug 2 ~u[ gi \( 4 - d y 2 , (28) 

X3 = u[gi} 4 

where 

g n ^[G{r)] n (29) 

r 

and where, for the usual reasons, g\ ~ (T — T c )~ l . 

This simple analysis already captures some qualitative 
features of the observed phenomena: a) It shows that 
generically the effect of even weak disorder grows stronger 
upon approach to criticality, and conversely that at tem¬ 
peratures well away from criticality the effects can be 
negligible as far as the initial growth of critical corre¬ 
lations is concerned, b) Moreover, it shows that gener¬ 
ically the first effect of disorder is to slow the growth 
of the susceptibility relative to its expected behavior in 
the absence of disorder. However, from this perturbative 
perspective, there is little reason to expect larger effects 
upon approach to a (T c = 0) quantum critical point than 
upon approach to a classical (T c > 0) critical point. 

Non-perturbative effects and quantum Griffith 
phases: However, it is well known that disorder - even 


weak disorder - can produce non-perturbative corrections 
to the thermodynamic properties of a system, and even 
change the analytic structure of the free energy by pro¬ 
ducing so-called Griffith singularities. For classical finite 
temperature systems, these Griffith singularities tend to 
be unobservably weak. However, for quantum systems 
with a discrete order parameter these effects are highly 
amplified producing a new set of phenomena known as 
quantum Griffith phases[50]. In particular, based on 
symmetry, we would tend to associate the universal crit¬ 
ical phenomena associated with a continuous quantum 
phase transition to a nematic ground state with those of 
the transverse field Ising model. 

The nature of the quantum Griffith phenomena associ¬ 
ated with the transverse field Ising model in the presence 
of random field disorder have been rigorously established 
in lDpU], but scaling arguments which are successful 
in accounting for these results produce compelling re¬ 
sults for higher dimensional problems as well[6lj- Among 
other things, these arguments imply that for a range of 
transverse fields h± greater than a critical value, h± >c , 
the susceptibility diverges in the limit T — > 0 as 

xCH ~ [1 /T} 1 - 6 with 0~ (hi- hi, c) (30) 

even though the ground-state does not exhibit nematic 
order. 

At core, the enormous difference between classical and 
quantum Griffith effects is illustrative of a more general 
tendency for the effects of disorder to be more dramatic 
on quantum critical phenomena than on their classical 
cousins. In classical statistical mechanics, a sum is per¬ 
formed over all position dependent field configurations 
with an appropriate Boltzman weight, and quenched dis¬ 
order implies the existence of point-like defects that lo¬ 
cally favor one particular value of the field. In quantum 
statistical mechanics, the same sum is performed over 
field configurations that are a function of position and 
imaginary time, and quenched disorder implies the ex¬ 
istence of line defects that are infinite in extent in the 
imaginary time direction. Thus, even rare impurities can 
have a large qualitative effect on the state of the system. 

It seems to us that the fact that the observed effects of 
disorder on the nematic susceptibility appear strongest 
where the extrapolated nematic ordering temperature is 
close to T=0 reflects the enhanced sensitivity of quantum 
critical phenomena to disorder. 

Non-equilibrium effects: Finally, it is important to 
note that non-equilibrium effects are unavoidable in any 
problem related to the random-field Ising model. The 
characteristic exponential slowing of dynamics upon ap¬ 
proach to the critical temperature implies that - neces¬ 
sarily - at a temperature strictly larger than T c , on any 
accessible laboratory time scale the system will fail to 
achieve thermal equilibrium. At all lower temperatures, 
glassy physics involving hysteresis, and noise must be 
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expected[12]. The existence or not of “aging” effects re- set of experiments, 
fleeting this physics has not been explored in the present 


Materials 

Fit 

range(K) 

2^66 

A/ao(K) 

T*(K) 

Sample 

dimension (pm) 

BaFe 2 As 2 Large 

135-250 

7.7±0.3 

-1540±13 

109.0±0.7 

3170x2920x50 

BaFe 2 As 2 Medium 

135-250 

7.3±0.8 

-1272±34 

116±2.2 

760x750x40 

BaFe 2 As 2 Small 

140-230 

5.4±0.3 

-948±13 

107.3±0.9 

280x300x20 

Ba(Feo.975Coo.025)2 AS2 

100-205 

14.5±0.8 

-2706±32 

77±0.8 

730x700x30 

Ba(Feo.953Coo. 047)2 As 2 

90-200 

9.4±0.4 

-2004±22 

59.9±0.9 

530x580x10 

Ba(Fe 0 .945CO0.055)2As 2 

70-180 

25.5±0.7 

-5370±44 

26.4±0.7 

580x540x25 

Ba(Feo.94Coo.o6)2As2 

50-197 

27.5±1.9 

-6812±158 

-1.4±2.7 

730x760x10 

Ba(Feo.93Coo.o7)2As2 

55-200 

13.7±0.5 

-3735±40 

1.8±1.1 

700x725x10 

Ba(Feo.92Coo.os)2As2 

50-200 

20.6±0.6 

-5758±60 

-26.1±1.4 

600x620x20 

Ba(Feo.9Coo.i)2As2 

38-200 

22.8±0.6 

-8773±99 

-94.4±2.1 

2950x2630x40 

BaFe2(Aso.68Po.32)2 

30-250 

7±1.9 

-1896F58 

11.7±3.1 

340x320x30 

BaFe2(Aso.64Po.36)2 

30-200 

6.6±1.3 

-2346±48 

-11.5±2.3 

280x300x10 

Ba(Feo.955Nio.o45)2As2 

60-210 

12.1±0.6 

-3487±44 

2.3±1.3 

750x700x15 

Bao.6Ko.4Fe2As2 

90-225 

-4.0±0.3 

820±19 

46.1±2.4 

750x760x25 

FeTeo.6Seo,4 

13-95 

-10.8±0.3 

-1400±17 

-11±0.7 

930x1100x30 


Table Sl:Fit parameters from the fit of 2m66 based on 2mee = 2mg 6 + \/[a(T — T*)] for all the compositions shown 
in Figures in the main text and supplementary material. 
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Figure SI: A schematic diagram illustrating the current and voltage configurations for measuring R xx and R yy using 
the modified Montgomery technique. A common ground point (A) is used for both configurations. The y Cartesian 
axis is aligned along the poling direction of the PZT stack. 
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Figure S2: A photograph and schematic diagram illustrating measurement of A p xx and A p yy elastoresistance for 
the i?2g symmetry channel of the parent BaFe2As2 compound. 





























18 




Figure S3:(a) A representative data set showing the measured change in resistance as a function of e xx strain. This 
data set was taken on optimally doped Ba(Fei_ a ;Co :E ) 2 As 2 (x = 0.07) at 240K. The zero e xx strain point is set at zero 
anisotropic strain, see text for details, (b) The normalized resistivity calculated from the raw resistance data shown 
in (a) using a modified Montgomery method. The normalization constant p is calculated at zero anisotropic strain as 
described by the main text. 
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Figure S4: (A p xx — Ap yy )/p as a function of anisotropic strain e xx — e yy at several temperatures for the optimally 
doped Ba(Fei_ a; Co a ;) 2 As 2 (x = 0.07). The normalization constant p is calculated at zero B 2 fl strain as described by 
the main text. Black lines show linear fits for each temperature, from which the elatoresistivity coefficient is 
extracted. Strains are positive for extension, and negative for compression. 
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Temp (K) 


Figure S5: Temperature dependence of strain transmission in the large BaFe 2 As 2 sample. A e yy and A e yy are the 
range of e yy strains experienced from a full voltage sweep applied to the PZT by the strain gauges mounted on 

Ae s 

the sample surface and PZT surface respectively. The ratio A v p is the fraction of strain transmitted through the 
sample, for perfect transmission this ratio is 1 and for no transmission the ratio is 0. Above 255K there is poor 
strain transmission which we tentatively attribute to glue softening and below 255K there is good (> 80%) strain 
transmission that shows a small temperature dependence. 
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Figure S6: A comparison between the -2m66 response of the large BaFe 2 As 2 sample calculated from the sample and 
PZT mounted strain gauge. The two curves are in good agreement with each other. Curie-Weiss fitting gives T* 
values of 109.OA' ± 0.7AT and 107. 7K ± 0.5A' respectively, indicating that the weak temperature dependence of the 
strain transmission below 255K does not affect the measured elastoresistivity coefficients. 
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Figure S7: Temperature dependence of the m.66 elastoresistivity coefficient for three different sizes of undoped 
BaFe 2 As 2 samples. The small sample (280/rm x 300/mi x 20/irn) has a smaller response than the medium (760/um x 
750/rm x 40/rrn) and large (3140/im x 3330/rm x 50/rm) samples. This implies equal strain transmission in the medium 
and large samples while a smaller strain is achieved for the small sample. 
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Figure S8: Comparison of the -2m66 elastoresistivity coefficient normalized by the value at 135A' for the large 
(3140/im x 3330/im x 50/im) and small (280/rm x 300^m x 20/im) samples. The two curves are in good agreement 
with each other and show a similar temperature dependence with T* values of 109 .OK ± 0.7 K and 107.7 K ± 0.5A' 
respectively. 
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Figure S9: Divergence of the i? 2 S elastoresistance 2777gg of several different families of optimally doped Fe pnictide and 
chalcogenide superconductors measured by the differential longitudinal elastoresistance technique using two separate 
crystals for each measurement. (A) optimally doped BaFe 2 (Aso. 7 Po. 3)2 (isovalent substitution), (B) optimally doped 
Ba(Feo.955Ni 0 .o45)2As2 (electron doped), (C) optimally doped Ba 0 .58K 0 .42Fe2As2 (hole doped), and (D) optimally 
doped FeTeo. 58 Seo. 42 - Insets indicate the dopant site (red) in the respective unit cells of each material. Upper 
panels show |2 to 66|, whereas lower panels show |(2 to 66 — 2mg 6 )|” 1 (left axes of lower panels, blue symbols) and 
|(2m 66 — 2m,g 6 )|(T — T*) (right axes of lower panels, black curves). Black(upper panels) and red(low panels) lines 
shows fits to Curie-Weiss behavior of 77755 and | (277155 — 277ig 6 )| _1 respectively. Grey horizontal lines (low panels) 
shows the average values of | (277755 — 277ig 6 )|(T — T*) in the fitting temperature range. Regions of deviation from 
Curie-Weiss behavior in (B) and (C) are indicated by gray shaded regions. For (A) and (B), 2?77gg is negative. For 
(C) and (D), 2mgg is positive. 
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Figure S10: Variation of the £? 2 9 elastoresistance of Ba(Fei_ a; Co a ;) 2 As 2 for four representative compositions measured 
by the differential longitudinal elastoresistance technique using two separate crystals for each measurement. (A,B) un¬ 
derdoped compositions Ba(Fe 0 . 975 Co 0 . 025 ) 2 As 2 and Ba(Feo.953Coo.o47)2As 2 ; (C) optimal-doped Ba(Feo.g 3 Coo.o 7 ) 2 As 2 ; 
and (D) overdoped Ba(Feo. 92 Coo.os) 2 As 2 . The heavily underdoped compositions is very well described by a Curie- 
Weiss temperature dependence over the entire temperature range(black lines in upper panels). For the compositions 
near optimally-doped, 2mee can be well fit by a Curie-Weiss T-dependence for temperatures at high temperatures. 
Below a characteristic temperature scale(different for each compositions, indicated by shade gray region), a strong 
downward deviation from Curie Weiss behavior is observed, also seen in the inverse susceptibility (upward curvature) 
and in |(2?7i66 — 2mg 6 )|(T — T*) oc Xn{T — T*) (strong downturn), which are shown in the lower panels. The deviation 
from Curie-Weiss behavior is the strongest at the optimal doping, and diminishes on either side of the phase diagram. 
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Figure Sll: Temperature-dependence of the resistivity for the same materials for which elastoresistance data are 
shown in the main text. The data were taken for free-standing (unstrained, and unattached) crystals. Data are 
plotted as R/R(300K) to eliminate uncertainty in geometric factors. 
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Figure S12: Temperature dependence of the elastoresistivity coefficients (mu — 77112 ) and 2m66 of (a) 
Ba(Fe 0 . 93 Coo.o 7 ) 2 As 2 (b) Bao.6Ko.4Fe2As2 and (c) FeTeo.6Seo.4. The data clearly show a divergence in the B 2g 
channel (i.e. 2m66 diverges) in all three iron based superconductors. In particular the elastoresistivity coefficients 
of FeTeo. 6 Seo .4 are similar to those of the Fe-arsenides, despite the fact that FeTe develops a spontaneous strain in 
the B\g symmetry channel at the structural phase transition. The data taken for this figure were from a differential 
elastoresistivity technique, see main text for details. 
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Figure S13: Temperature dependence of the normalized deviation 6 of A/a for Ba(Feo.o 93 Coo.o 7 ) 2 As 2 . <5 is calculated 
by taking the difference of (2?ri66 — 2m/ 6 )(T — T*) from its mean and then normalized by the mean based on the fitted 
2 ?7ig 6 and T*. The fitting ranges from top to bottom are: 25 - 300K, 25 - 160K, 55 - 160K, 55 - 205K. Arrow in the 
top panel indicated the deflection point of high temperature deviations. 
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Figure S14: Temperature dependence of the 1?2 S elastoresistance coefficient 2m66 of Ba(Fei_ a; Coi C ) 2 As 2 for eight 
doping concentrations from a: = 0 to a; = 0.1. Solid curves shows fits to the Curie-Weiss behavior, and the fitting 
parameters are tabulated in table SI. Data were successively offset by 10 for clarity. 








30 



(T - T*) (K) 


—P optimal 
— Ni optimal 
—Co optimal 
—■ K optimal 
—Fe(Te,Se) optimal 


Figure S15: Power law behavior of the B 2 g elastoresistance of optimally doped BaFe 2 (Aso. 7 Po. 3 ) 2 - Figure shows 
log[|(2TO66 — 2wg 6 )|] v.s. log (T — T*) for all measured optimally-doped iron-based superconductors. 2mg 6 and T* were 
extracted from the Curie-Weiss fit [36]. Only the optimally P-substituted data (red) can be fitted linearly over the 
whole temperature range, with a slope 7 = 0.985 ± 0.005 (solid black line). Note that for Co, Ni, and P substitution, 
2 to 66 ^ 2?rig 6 is negative, whereas for K substitution and FeTei-^Sex, it is positive. 
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Figure S16: Temperature dependence of the £? 2 S elastoresistance coefficient (a) 2m66 and its inverse (b) |(2 to 66 — 
of BaFe 2 (Asi_ a ;P a ;)2 for x = 0.32 and x = 0.36. Solid curves shows fits to the Curie-Weiss behavior, and 
the fitting parameters are tabulated in table SI. From, the x-intercepts of the linear fits of |(2m66 — 2TOg 6 )| _1 it can 
be seen that the Weiss temperature crosses zero as the doping concentration increases from optimal doped to slightly 
overdoped. 









